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Abstract. Results from direct detection experiments are typically interpreted by employing 
an assumption about the dark matter velocity distribution, with results presented in the 
m-^ — Un plane. Recently methods which are independent of the DM halo velocity distribution 
have been developed which present results in the Vmin — 9 plane, but these in turn require an 
assumption on the dark matter mass. Here we present an extension of these halo-independent 
methods for dark matter direct detection which does not require a fiducial choice of the 
dark matter mass. With a change of variables from Vmin to nuclear recoil momentum (pr), 
the full halo-independent content of an experimental result for any dark matter mass can 
be condensed into a single plot as a function of a new halo integral variable, which we 
call h{pR). The entire family of conventional halo-independent g{vmin) plots for all DM 
masses are directly found from the single h{pR) plot through a simple rescaling of axes. By 
considering results in h{pR) space, one can determine if two experiments are inconsistent 
for all masses and all physically possible halos, or for what range of dark matter masses 
the results are inconsistent for all halos, without the necessity of multiple g{vmin) plots for 
different DM masses. We conduct a sample analysis comparing the CDMS II Si events 
to the null results from LUX, XENONIO, and SuperCDMS using our method and discuss 
how the mass-independent limits can be strengthened by imposing the physically reasonable 
requirement of a hnite halo escape velocity. 


Contents 


1 Introduction 1 

2 Review of Halo-Independent Methods 3 

3 Halo-Independent Direct Detection Without Mass Assumptions 6 

3.1 Direct Detection in Momentum Space 6 

3.2 Comparisons Between Null and Positive Results 8 

3.3 Additional Constraints from Escape Velocity 12 

4 A Sample Analysis: Comparing CDMS H Si with SuperCDMS, LUX, and 

XENONIO 12 

5 Analysis Recipe 15 

6 Conclusions 16 

A Comparisons Between Null Results 17 

A.l Ratio Plots for Null Results 17 

A.2 A Sample Analysis: Comparing LUX, XENONIO, and SuperCDMS 19 


1 Introduction 

The strong astrophysical evidence for dark matter (DM) has prompted an extensive program 
of experiments to detect DM terrestrially [1-9]. Such direct-detection experiments exploit 
the possibility that DM can scatter off nuclei in a large detector volume [10], with a DM 
scattering event identified through the ensuing detection of the nuclear recoil energy Er. 
However, the event rate as a function of is strongly dependent on the incident DM 
velocity distribution. In particular, the scattering rate is proportional to the integral over 
the entire velocity distribution of dark matter in the Galactic halo, which has never been 
measured directly. 

Eor simplicity it is often assumed that the speed distribution is Maxwell-Boltzmann in 
nature. However, even in this simple case varying the Sun’s circular velocity and the Galactic 
escape speed can significantly alter the interpretation of direct detection results [11-13]. Sim¬ 
ilarly, marginalising over the parameters associated with the speed distribution motivated by 
dark matter numerical simulations, which show departures from pure Maxwell-Boltzmann, 
leads to uncertainty in interpretation of results [14, 15]. Given the coarse resolution of numer¬ 
ical simulations and the uncertainties associated with the extracted speed distributions, one 
may instead parameterise the speed distribution in some general way [16-19]. With sufficient 
data, and a judicious choice of basis for the parameterisation, one can hope to determine both 
the astrophysical and particle physics properties of dark matter simultaneously [18, 20]. 

However, motivated by our present lack of understanding of the local dark matter ve¬ 
locity distribution, techniques have been developed which interpret direct detection results 
independent of assumptions about the astrophysics. Given at least two experimental mea¬ 
surements of dark matter scattering, on different elements, it is possible to determine the dark 
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matter mass without making assumptions about the form of the velocity distribution [21], 
Even without two positive results it is possible to compare experiments in an astrophysics 
independent fashion by constraining the halo integral itself, rather than the DM-nucleus scat¬ 
tering cross section [22]. For a given nuclear recoil energy, only DM traveling at sufficiently 
high velocities can provoke a nuclear recoil at that energy. Thus, such halo-independent con¬ 
straints are a function of the minimum DM velocity VminiER) required to provoke a recoil 
energy Er, which for elastic scattering is 

VminiER) = . (1.1) 

While we will limit our discussion to the case of elastically scattering dark matter in a single 
target experiment, these astrophysics-independent techniques have been extended to include 
multiple targets [23], inelastic dark matter [24], momentum-dependent scattering [25], other 
more general forms of dark matter scattering [26], searches for the modulating signal [27], 
and have been applied to compare data from many experiments [28-36]. The techniques 
have also been extended to enable a halo-independent combined-likelihood function to be 
constructed [37] , allowing for the tension between DM hints and null results to be quantified 
in a systematic way, see also [38]. 

As can be seen from Eq. (1.1), Vmin depends explicitly on the DM mass through 
the DM-nucleus reduced mass fiNx- Thus, to analyze a direct detection experiment using 
existing halo-independent methods, one must make a choice for rriy-. This is undesirable 
since it leads to a proliferation of plots, one for each choice of mass. Furthermore, the goal 
of halo-independent methods is to “factor out” uncertainties about dark matter properties. 
Indeed, the DM mass is probably less well-constrained than the velocity distribution: while 
DM velocity in the Galactic halo is bounded by the Galactic escape velocity, there exist 
well-motivated DM candidates with masses spanning some twenty orders of magnitude [39]. 
Even in the range of masses 1 GeV — 1 TeV typically probed by direct detection experiments, 
different choices of DM mass can alter constraints derived from halo-independent analyses, 
and hints of positive signals can be either excluded or allowed in different regions of Vmin 
depending on the DM mass [35]. 

In this paper we show that, for DM which undergoes non-relativistic elastic scattering 
with nuclei, many ambiguities related to the unknown DM mass can be resolved by a simple 
change of variables. By presenting limits and preferred values for the halo integral as a 
function of nuclear recoil momentum pR, instead of DM minimum velocity Vmin or nuclear 
recoil energy Er, one can condense the halo-independent results from an experiment in a 
way which is valid for any value of the DM mass. This is a direct consequence of a kinematic 
relation for non-relativistic elastic scattering, derived from Eq. (1.1), 

PR = ‘^PNxVminiER,) = \/ 2mArE/j , (1.2) 

which implies that pR is related to Vmin by a constant rescaling, but is manifestly independent 
of To put it another way, relating the observable Er to Vmin requires a choice for the 

^This is the essential reason for working in pjj-space. A similar statement cannot be made about an Er- 
space plot where, for nonrelativistic scattering, the rescaling of the x-axes into Vmin-spa.ce is nonlinear in Er 
and explicitly depends on the DM mass. For relativistic scattering, though, one can make a mass-independent 
plot in Tij-space [40]. 
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DM mass, whereas relating Er to pR does not.^ Thus, for a given experiment, the entire 
1-parameter family of Vmin plots for all DM masses can be condensed into a single pR plot, 
where the scale on the pR-axis does not change as is varied.^ Furthermore, pR enjoys the 
same advantage of Vmin in that multiple experiments with different nuclear targets, which 
may have different ranges of sensitivity in Er, can in principle have overlapping sensitivity 
in PR. 

While momentum-space plots are useful for a single experiment, the dependence of pR 
on the target mass rriN means that care must be taken in comparing two different experiments 
on the same pR plot. Naively one would expect that the comparison must be carried out on a 
mass-by-mass basis, necessitating again that a large number of plots be made and compared. 
However, we show that this issue can be avoided entirely by considering the ratios of exclusion 
curves or envelopes of preferred values in p/j-space. By plotting these ratios as a function of 
a rescaled halo integral h, which does not depend on the target mass, a single plot contains 
sufficient information to determine if a dark matter signal is excluded for all dark matter 
halos and for all dark matter masses, or whether there is a range of masses for which there 
is agreement. The rescaled halo integral h has a conceptually straightforward interpretation 
as the differential scattering rate normalized to a single nucleon. In the most general case, 
it is simple to read off from this plot the range of DM masses for which two experiments are 
inconsistent for any DM halo. 

This paper is organized as follows. In Sec. 2, we review the general setup of halo- 
independent methods, with particular emphasis on the role of Vmin- In Sec. 3, we describe 
the change of variables from Umm-space to pR-space and show how this can be used to present 
experimental data in a halo-independent manner which is also independent of the DM mass. 
The key results of this paper are in Sec. 3.2, where we show how p/j-space can be used to 
compare null results from one experiment with positive results from another experiment to 
derive regions of inconsistency which are valid for all DM masses above or below (depending 
on the ratio of target masses) a critical value. Finally, in Sec. 4, we perform a sample analysis 
on CDMS II Si, LUX, XENONIO, and SuperCDMS data to illustrate the utility of p^j-space. 
For the convenience of the reader we summarize the minimal analysis recipe in Sec. 5. We 
conclude in Sec. 6. In App. A we comment on the use of mass-independent methods for 
comparing the relative strengths of different null results. 

2 Review of Halo-Independent Methods 

Here we briefly review the standard techniques for halo-independent analyses of direct de¬ 
tection experiments. For spin-independent DM-nuclear scattering, the differential event rate 
at a direct detection experiment is given by 

^IEr ~ 2m jjfl J dE'j^G{ER,E'pi)e{E'pi)E‘^{E'pi)g{vmin{E'ji)) , (2.1) 

where is the DM mass, the nucleon mass, pnx nucleon-DM reduced mass, an 
the DM-nucleon scattering cross-section, the local DM density, Na is Avogadro’s number, 

^In [18] it was also previously emphasized that the recoil momentum is an attractive variable for analyzing 
direct detection data, especially as it is independent of the DM mass. 

^This conclusion holds only for elastic scattering. The extent to which one can derive constraints without 
mass assumptions for more general kinematics is explored in a forthcoming paper [41]. 
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F(Efi) is the nuclear form factor which accounts for loss of coherence as the DM resolves sub- 
nuclear distance scales, Ct{A,Z) = (fp/fnZ -|- (^ — Z)) is the usual coherent DM-nucleus 
coupling factor, e{Eji) is the detector efficiency, and G{Er, E'j^) is the detector resolution 
function. 

The integral over the velocity distribution (which we will often refer to as the “halo 
integral”) is 





V 


( 2 . 2 ) 


where f{v) is the unknown DM velocity distribution and ve the Earth’s velocity, both 
in the Galactic frame. We shall ignore the small time dependence of the Earth’s velocity in 
the Galactic frame. The lower limit of the halo integral, Vmin, is the recoil energy-dependent 
minimum DM velocity required to produce a nuclear recoil and depends on the kinematics 
of the interaction. Eor elastic scattering Vmin{Eji) is given by Eq. (1.1). 

As first pointed out in [22, 42], since the DM velocity distribution f{v) is positive semi- 
dehnite, g{vmin) must be a monotonically decreasing function of Vmin- This leads to powerful 
constraints on the shape of the halo integral derived from either null or positive experimental 
results, which we will review shortly. Eirst, though, we summarize the essential differences 
between the usual — an analyses and halo-independent analyses. 


• m — a plots. Given a choice of DM velocity distribution f{v) and DM local density 
event rates can be calculated as a function of and an from Eq. (2.1). Null results 
give exclusion contours, while positive results generally give closed preferred regions in 
— an space. Eor each choice of f{v), a different plot must be made. 


• Halo-independent plots. By rescaling the halo integral as 


di^min) 


Px'^n 




(2.3) 


to absorb all the detector independent terms in Eq. (2.1), event rates can be calculated 
from Eq. (2.1) as a function of g, making the most conservative choice for f{v). A 
choice for is still necessary to relate Vmin to Er through Eq. (1.1). Null results 
once again give exclusion contours, but in Vmin — g{vmin) space. Positive results give 
preferred shapes for the halo integral g{vmin), which is constrained to be a monoton¬ 
ically decreasing function; thus positive results generally give open envelopes of such 
curves in Vmin — g{vmin) space, rather than closed regions. Eor each choice of a 
different plot must be made. 


The choice of Vmin as the x-axis variable for halo-independent plots is one of convenience, as 
it is the natural argument for the halo integral g, but since VminiEn) is monotonic for elastic 
scattering, one may just as easily present results in Er — p(Er) space, since g{EE) is also 
monotonically decreasing. The advantage of Vmin is that experiments with non-overlapping 
sensitivity in Er may have overlapping sensitivity in Vmin-, due to the dependence of Vmin 
on the target mass mjv- However, we will see in Sec. 3 that the choice of pR for the x-axis 
variable shares this same advantage, while at the same time obviating a choice of m^. 

To derive halo-independent exclusion curves, note that for a given point vq on the Vmin- 
axis, the halo integral which gives the least number of events consistent with monotonicity 
is 

g(Fmin) — 5O0(^O Vmin) 
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Setting an upper bound on the number of events using this choice of halo integral, one derives 
the weakest {i.e. most conservative) possible bounds on g{vo)- Thus, by sweeping over vq in 
the Vmin range of an experiment, one can build up an exclusion contour which at every x- 
coordinate vq sets the most conservative bounds on the halo integral: the true halo integral 
must lie entirely below this exclusion contour for every Vmin- Such an exclusion curve is 
extremely robust, as any positive signal so excluded can be said to be excluded for all DM 
halos. 

Determining best-fit values and envelopes for the halo integral for a positive result 
can be done in many ways. As a representative example, one may interpret Eq. (2.1) as a 
component of an unbinned extended likelihood function 


C 


^ dR 
No\ l^dER 


En—Ei 


(2.5) 


where Nq total events are seen at recoil energies The monotonicity of g can then be 
used to reduce the possible shapes of g{vmin) to a finite-dimensional subset, allowing a simple 
numerical maximization of the likelihood function. For details, we refer the reader to [43] 
where these statements are proven and expanded upon, as well as [44] which describes an 
alternate maximization technique. Including finite energy resolution effects, the preferred 
region at a given confidence level is an envelope in Vmin — g{vmin) space. Thus, for any point 
inside the envelope there is a curve which lies entirely within the envelope and which passes 
through that point which is consistent with the data at the required confidence. In addition, 
there are no curves within the required confidence containing points which lie outside the 
envelope. Consequently, an exclusion limit which crosses the lower boundary of this envelope 
is able to rule out all DM halos preferred by the data, and thus excludes the putative signal 
in a halo-independent fashion. 

When carrying out a halo-independent analysis, as mentioned above, a particular DM 
mass must be assumed. The constraints on, or predictions for, g{v) for another DM mass can 
be determined by simple rescalings [43]. For a DM mass, a point {vmin,g) is mapped to 
a new point {v'^i^,g') for DM mass by 

{Vmin 1 g 1 ^x) ^ f '^min ; n 5 ) ^x' 1 ' (^•®) 

\mx' l^nx j 

Since the mapping depends on the mass of the nuclear target, this will shift different detectors 
by differing amounts. In the following section we will present a change of variables which 
allows the presentation of results in a fashion independent of the DM mass. Essentially this 
change of variables ‘factors out’ the re-scaling freedom of Eq. (2.6) to leave a single plot 
which is the same for all m^. 

Of course, one could alternatively present results in the usual Vmin — g{vmin) space plot 
for a single choice of rriy^ and then all other Vmin — g{vmin) space results for different masses 
could be immediately found from the rescaling of Eq. (2.6).® However, there would be am¬ 
biguity surrounding which fiducial DM mass to take when presenting results in publications. 
Furthermore, as we will see, the pR — h{pR) space option advocated here allows for the po¬ 
tential inconsistency of two different experimental results for any halo and any DM mass to 

^For large No the events may be binned. In this case it is very simple to extract the preferred region of 
g{vmin) from binned likelihood or methods. 

®This was described in more detail in [43]. 
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be determined directly by eye from the pR — h{pji) space plot, whereas this is not possible 
for a Vmin — g{vmin) space plot with a specific choice of DM mass. 

3 Halo-Independent Direct Detection Withont Mass Assnmptions 

In Sec. 3.1 we will describe the change of variables for calculating momentum-space con¬ 
straints for a given experiment and then discuss how the standard halo-independent plot for 
any choice of DM mass may be easily constructed from this single plot. The momentum- 
space plot is unique as only one such plot exists for any set of direct detection data, and is 
universal in the sense that it applies for all dark matter masses, and does not need to be 
reconstructed on a mass-by-mass basis, as is the case for the usual Umm-space plots. Follow¬ 
ing this, in Sec. 3.2 we demonstrate the utility of the momentum-space plot by showing how 
the momentum-space constraints may in some cases be used to draw strong and general con¬ 
clusions on the comparison between different direct detection experiments which hold for all 
DM halos and for all DM masses, making these conclusions halo-independent and true for all 
masses. Alternatively, the same information may, in some cases, imply a range of DM masses 
for which experiments are inconsistent, independent of DM halo. We show in Sec. 3.3 how 
these conclusions may be strengthened by restricting to halos with a finite escape velocity. 

3.1 Direct Detection in Momentum Space 

Eq. (2.1), which describes the differential DM event rate, contains a number of unknown fac¬ 
tors which depend on the DM mass, cross section, and local density. To compare observations 
with specific models, one must make choices for these factors, but to simply compare results 
at different detectors this is unnecessary. Thus, we may absorb all detector-independent 
factors into a single unknown, which we write as 

HPr{Er)) = g{Vmin{ER)) = g{Vmin) ■ ( 3 . 1 ) 

We emphasize that the rescaling factor between h and g is independent of the detector 
nucleus, since it only depends on the nucleon mass nin- As our notation makes clear, we are 
thinking of h as a function of recoil momentum, 

PR = \J 2mNER = 2pMxVmin , (3.2) 

which for a given recoil energy Er is independent of DM mass, rather than Vmin which is 
not. Both h and g are implicitly monotonically-decreasing functions of Er, since Vmin is a 
monotonic function of Er and the halo integral g is constrained by the arguments of Sec. 2 
to be a monotonic function of Vmin- We can now rewrite Eq. (2.1) as 

= CUA,Z) I dE'^G{ER,E'R)e{E'R)E\E'^)h{pR{E'^)) . (3.3) 

Notice that in Eq. (3.1) we break with the convention of [22, 42], and instead of the quantity 
g we opt for h which absorbs all dependence on the DM mass, and is also independent of 
the detector. The use of h also has a conceptual advantage because it is a physical quantity: 
the differential scattering rate normalized to a single nucleon, before correcting for resolution 
and detection efficiency, analogous to the normalized WIMP-nucleon cross section used by 
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experimental collaborations to report limits. This interpretation makes clear that halo- 
independent comparisons of experiments are just direct comparisons of recoil energy spectra, 
after a simple rescaling to account for the kinematic effects of different target masses.® 

Since there is a one-to-one relation between Vmin and it is simple to repeat the 
discussion of Sec. 2 in terms of pR. For instance, if DM scattering were to produce an event 
with nuclear recoil momentum po j then the halo integral which could lead to this event while 
producing the minimum number of scattering events at other recoil momenta is described by 

Hpr) = hoQiPo - Pr) ■ (3.4) 

Here, po is taken over the whole range of pR sensitivity of the experiment: if the experiment 
is sensitive to energies satisfying Eiow < E < Eugh, then ^/2misfEiow < PR < 

Thus, one can set limits on /iq, and thus /i(po), based on the null results of an experiment by 
integrating over Er and e.g. using a Poisson upper limit on the total number of events at a 
given confidence level, for each value of po. Just as easily, h{pR) can be chosen to be a sum 
of step functions as proposed in [43, 44] in order to find best-fit regions for the halo integral 
in the case of positive signals. In either case, by construction this analysis makes only the 
most conservative assumptions about the DM velocity distribution. We emphasize that at 
no point in calculating the constraints on Eq. (3.3) is it necessary to make a choice for m^. 

It is convenient at this stage to make a few observations on the possible form of 
momentum-space exclusion curves which will be useful later. 

• As a detector is not sensitive to nuclear recoils at vanishing momentum, and most 
analyses have a low-momentum threshold, the exclusion curves will in general become 
weaker at low recoil momenta and asymptote to infinity at some fixed small recoil 
momentum. 

• As the most conservative choice of halo (3.4) predicts larger numbers of events for 
larger recoil momenta, the limits will in general become stronger towards higher recoil 
momenta. The result is that the exclusion curve in momentum-space will in general be 
a monotonically decreasing function. ‘ 

• Detectors in principle have sensitivity out to infinite recoil momenta, as a signal beyond 
the high momentum threshold of a detector implies a nonzero number of predicted 
events within the signal region. Thus, the momentum space exclusion contour will in 
principle extend to infinite recoil momenta, but due to the finite energy range of the 
signal region, the bound will approach a constant value in h. 

The momentum-space representation of experimental results contains all of the infor¬ 
mation conveyed by the usual Umm-space plots, and these plots may actually be constructed 
directly from one another. From Eqs. (3.1) and (3.2), the mapping is linear for any choice of 
DM mass: 

(“’ *) ^ 

®We thank Felix Kahlhoefer for emphasizing this interpretation to us. 

^Although h should, on general grounds, be a monotonically decreasing function, the limits on h extracted 
from real data may not be monotonically decreasing. For example, if Yellin’s methods [45] are used, then 
limits can increase near the Vmin where the limit-setting method switches from one interval to another. 
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Figure 1. (a) A schematic of the momentum-space plot proposed in this work for two different 
detectors with target nuclei masses satisfying Ma > Mb- (b) A schematic of the standard Umin-space 
plot. The mapping from the unique momentum-space plot to the Umm-space plot is also shown in (a) 
for a typical choice of DM mass making clear that the differences only arise in the horizontal direction, 
see Eq. (3.5). 


From the behaviour of g{v) under changes of DM mass, Eq. (2.6), we see that the single 
momentum-space plot contains all halo-independent information. 

We can also compare the results of two different experiments on the same pji plot. For 
any DM mass, the mapping of Eq. (3.5) rescales both curves in the vertical direction by the 
same factor, but rescales them in the horizontal direction by different factors depending on 
both the DM mass and the target nucleus mass. Hence, in mapping from recoil momentum 
space to Umin-space, exclusion curves for two different detectors are shifted relative to each 
other only in the horizontal direction, as depicted in Fig. 1. We now show how this may 
be exploited to draw very general halo-independent conclusions on the relative behavior of 
different experiments which hold for all DM masses. 

3.2 Comparisons Between Null and Positive Results 

A principal purpose of the momentum-space method will be to serve as a general tool to 
determine the consistency of potential positive results with bounds from detectors which 
observe no hint of DM scattering. There is a great deal of value in this sort of comparison, as 
it may shed light on the nature of an emerging DM discovery and possibly point towards DM 
mass ranges which are preferred for such a signal. In general, positive results give preferred 
regions which explain the observed scattering events in a detector to within some given 
statistical confidence. Technology for determining these best-fit regions has been developed in 
the context of halo-independent Umm-space plots and may be easily adapted for momentum- 
space plots. As an example, using the method described in [43] it is possible to determine 
the preferred parameter space using an unbinned likelihood function, leading to preferred 
regions such as those depicted schematically in Fig. 2. Unlike the preferred regions in m^ — an 
space, the preferred parameter regions in Umm-space will not in general be closed due to the 
monotonicity constraint on the halo integral. 

At this point it is worth discussing some general features of the preferred parameter 
regions in momentum space. The upper boundary of a preferred region will be very similar 
to an exclusion contour. Signal events are preferred, but not too many; this boundary acts as 
an exclusion contour which may extend to non-zero values for arbitrarily high recoil momenta. 
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Figure 2. Illustrative possibilities of momentum-space plots for results from an experiment B (shaded 
red) which observes a tentative DM signal and a null experiment A (black lines) which is used to set 
constraints. The target nuclei masses are here assumed to satisfy Ma > Mb- The rescaling to map 
from momentum-space to Wmm-space is only shown in one case for clarity, but this mapping holds for 
all cases. In (a) an experiment A finds constraints which exclude a DM interpretation of all anomalous 
events observed in B as the exclusion contour crosses the lower boundary of the preferred region. Any 
mapping to Umin-space will shift the region for B further to the right relative to A, implying that 
experiment A excludes a DM interpretation of B for all DM masses and all DM halos. In (b) it may 
be that A excludes a DM interpretation of B for all halos for some finite range of DM masses. In (c) 
and (d), experiment A can never exclude a DM interpretation of B for all halos for any choice of DM 
mass, as shifting B to the right never causes A to cross the lower boundary of B. However, in (c) 
there is tension between the results of A and B independent of the DM mass, whereas in (d) there 
is only tension for a range of DM masses. See the text for a general method to resolve and quantify 
these differences. In green we also show the construction of the variable R{h), defined in Eq. (3.7) 
which allows for a quantitative assessment of all possibilities. 


It also becomes weak below the threshold of the detector, approaching infinity as the recoil 
momentum tends to zero. The lower boundary of a preferred region is different. Above the 
high energy threshold of a detector there is no preference for signal events, thus the lower 
boundary will tend towards zero at high recoil momenta. At recoil momenta below the low 
energy threshold there is also no preference for signal, however monotonicity combined with 
the preference for signal at higher momenta means that the lower boundary of the preferred 
region will approach a fixed value at low recoil momenta, which we denote hiow All of these 
features are depicted in Fig. 2; it is worth emphasizing that these features would be expected 
for a momentum-space plot for any DM hint. 

We consider two experiments: the null experiment A which may be used to place con¬ 
straints in momentum-space, and experiment B which observes anomalous scattering events 
which may be interpreted as DM scattering. In the schematic of Fig. 2, the black curves 
are from A and the red shaded regions are from B. In some simple cases it is possible to 
determine the relative strengths of bounds from two different detectors directly from their 
hipn) plots as described in the caption. 

In more general cases, determining the outcome of mapping the momentum-space exclu¬ 
sions to Umm-space for a specific DM mass is less straightforward. To treat the general case we 
exploit the fact that mapping from momentum-space to Umm-space only shifts curves relative 
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to each other in the horizontal direction. Each momentum-space exclusion curve or lower 
boundary of a preferred region is given by a function h{pji); since this curve is monotonically 
decreasing, we may define the inverse mapping PR{h) over the fixed range 


h £ [h{pR —>• oo), oo] 


(3.6) 


In order to quantify the relationship between exclusion contours and preferred regions of 
parameter space it is useful to consider ratios of pR as a function of h. Specifically, we define 


R{h) = 


PR,low (^) 


(3.7) 


where PR^nm is an exclusion curve, and Pr,iow(^) is the curve describing the lower boundary 
of a preferred region.® This construction is shown in green in Fig. 2. In the case of a positive 
signal, it is the lower boundary of the preferred region which determines the consistency 
between a DM interpretation of experiment B relative to an exclusion from experiment A. 
We also define a related quantity 


R^{h) = 


^min,lim(^) 


'^min,\aw 


Ch) ’ 


(3.8) 


which is the ratio of the two curves in Vmin space.® As the notation makes clear, R^ depends 
explicitly on m^. If Rx{h) < 1, the exclusion contour from A crosses below the lower 
boundary of the preferred region from B, and we can say that the signal from B is inconsistent 
with the null results from A for the particular mass which defines the mapping into Vmin- 
space. 

The relationship between R and R^ is a simple scaling. 




with 


Fx = 


PBx _ Ma + Mx 


(3.9) 


(3.10) 


PAx Mr + Mx 

This scaling is inherited from the scaling of pr given in Eq. (3.5). Crucially, over all DM 


masses satisfies the inequality 


^<Fx<l, A Ma> Mb 


1<F 


X < if Mb > Ma 


(3.11) 


This restricted range of Fx is a key tool for extracting mass-independent statements concern¬ 
ing different experiments. Consider the case of Ma > Mr- If R lies entirely above Ma/Mr, 
then the two results are consistent for all DM masses, since as the DM mass is varied, can 
never dip below 1 thanks to the inequality (3.11). If any of the R curve lies below 1, then A 
and B are inconsistent for all DM masses since some point on Rx will lie below 1 for all m^. 
Similar reasoning holds for the case Ma < Mr, but with the roles of the critical lines at 1 
and Ma/Mr reversed; see Fig. 3. 


®In some cases, such as (c) and (d) of Fig. 2, R is not well-defined as there is no overlap in the regions of 
h explored by the exclusion contour and the lower boundary of the preferred region. In these cases it is clear 
by eye that there are no DM masses for which A can exclude B independent of the DM halo. 

® Since the mapping from h to g is the same for both A and B, we prefer to think of as a function of h 
rather than g for simplicity of notation. 
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Ma>Mb 


Mb>Ma 




h h 

Figure 3. Some possibilities for the variable Eq. (3.7), for both possible mass orderings. When 
moving from p/j-space to r;mm-space, the height of any line is rescaled by the DM mass-dependent 
parameter Eq. (3.10) which satishes the appropriate inequality in Eq. (3.11). Case (a) in Fig. 2 
corresponds to the solid red curve, and case (b) corresponds to the dashed or dot-dashed red curves. 
Cases (c) and (d) cannot be plotted because the lower boundary and the exclusion curve have no 
overlap in h, and R{h) is undefined. 


Suppose the global minimum of R{h) is located at h = Hq. In the case where R{ho) lies 
between 1 and Ma/Mb (for either mass ordering), define 


5 = R{ho) - 1 


(3.12) 


Note that 5 < 0 for the case Mb > Ma- Depending on S, there will be a range of DM masses 
for which A excludes B. The lower or upper endpoint of that range (depending on the mass 
ordering) will be the value of where i?^(/io) = 1, given by 


Mab{S) 


MaMbS 
Ma-{1 + 5)Mb 


(3.13) 


For the two cases the range of DM masses for which A excludes B is 


r > Mab(< 5), if > Mb . . 

\m^< Mab{S), if Mb>Ma ■ ^ ’ 

A schematic of these various scenarios is shown in Fig. 3. 

We emphasize that all of the results discussed above hold for all DM halos because the 
constraints are constructed for the most conservative halo possible. Furthermore, carrying 
out the mapping into VminSpace to obtain R^ is actually unnecessary, since the R{h) plot 
contains all the necessary information. Thus in practice, the entire analysis can be carried 
out without a choice of by working exclusively with R{h)}^ 

We focus here on pairwise comparisons of one hint (defined by a best fit region at 
confidence level CLhint) with one exclusion (dehned by a bound at confidence level CLbound) 

^°The mapping into Dmin-space in this section was performed to make explicit the connection between 
previous methods and this method. 
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at a time. It would be straightforward to combine all hints in a joint likelihood function and 
compare this region with each constraint individually. We do not attempt to quantify the 
level of tension, if there is any, beyond the statement that the two results are inconsistent at 
their respective confidence levels. Such techniques do exist [38], but we have not yet tried 
to extend them to /i-space. We have assumed the bounds are derived using techniques that 
do not require the form of the background be known, such as the techniques of Yellin [45] . 
However, if the background is well-modeled one can go further and combine both signals 
and bounds into one global likelihood function and compare the relative likelihood of various 
hypotheses to fit all the data [37] . 


3.3 Additional Constraints from Escape Velocity 


Up to this point, our analysis has been completely halo-independent, in the sense that we have 
imposed no constraints on the DM velocity distribution, no matter how physically reasonable. 
Astrophysical halos, however, have escape velocities, and it is worth considering the effect 
of imposing a cutoff on the velocity distribution at Vesc on the behavior of momentum-space 
exclusion limits. In Umm-space, a cutoff at Vesc means that a preferred region for a positive 
signal at Vmin > Vesc should no longer be interpreted as coming from DM scattering. Consider 
again the relation between Vmin and pR, 




PR 

2flN), 


(3.15) 


Thanks to the reduced mass in the denominator, the rescaling to Umin-space is a monotoni- 
cally decreasing function of m^, which means that for sufficiently small m^, a given point on 
the pR-axis will always be mapped to Vesc- 




rriNPR 

2mNVesc — PR 


PR ^ Vesc 


(3.16) 


For this pR value, only larger than can be consistent with a DM interpretation. 

As a consequence, halo-independent exclusion limits which hold only for limited DM 
mass ranges may in fact hold over all DM masses consistent with physically reasonable halos. 
Suppose that Ma > Mr, and let PB,max be the largest pij-coordinate of the lower boundary 
of envelope B. Setting rriM = Mr in Eq. (3.16) and comparing with Eq. (3.14), one can 
conclude that A excludes B for all DM halos with escape velocity less than or equal to 


{Ma - Mr)pr 

2MaMr6 


(3.17) 


Put another way, the minimum DM mass required to bring A and B into agreement is already 
small enough that every halo in the preferred region implies Vmin values greater than this Vesc- 
Note that such an analysis is not possible when Mr > Ma, because in that case > 1 and 
the relevant inequalities are reversed. However, if m^^min > Mar (6), then the constraints 
on DM masses from the escape velocity are stronger than the requirements for consistency 
between results from R{h). 


4 A Sample Analysis: Comparing CDMS II Si with SnperCDMS, LUX, 
and XENONIO 

To demonstrate the methods of section 3.2 for comparing positive signals with exclusions, 
we consider the case of the SuperCDMS, LUX, and XENONIO limits and the preferred halo 
for the Si detectors of the CDMS H experiment. 
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The SuperCDMS collaboration has publicly released data from their analysis of 577 kg d 
of exposure on their seven Ge detectors with the lowest energy thresholds [7, 46] . The blind 
analysis observed 11 events with GT^'^Qg events expected from background. Although the 
observed events are consistent with the background estimate, the SuperCDMS collaboration 
did not perform a background subtraction before calculating the limit, conservatively as¬ 
suming that all observed events are potential DM scattering. Using the efficiency, ionization 
yield, and energies of the 11 events from the data release, we set 90% C.L. limits using the 
“Pmax” method [45]. 

The CDMS II Si analysis observed three events in 140.2 kg d of exposure. We follow the 
analysis approach used in [43], taking a gaussian detector resolution of 0.2 keV, the acceptance 
from [47], and backgrounds from [48]. We furthermore assume that the distribution of AL is 
corresponding to five degrees of freedom: one for each step position and height in the halo, 
minus one for the monotonicity constraint. The assumption is a reasonable approximation 
to the sampling distribution of AL [43], a more precise determination could be made by a 
Monte Carlo simulation, varying the number of events in each pseudoexperiment. A value of 
AL = 9.2 then corresponds to 90% C.L. We perform Markov Chain Monte Carlo (MCMC) 
sampling^^ of the 6-dimensional parameter space, and use halos with AL < 9.2 to calculate 
the envelope preferred by the CDMS II Si result. 

For comparison, we also compute limits for LUX and XENONIO at a fixed dark matter 
mass as described in section 2. Similar to [43], we use events and ionization yield Qy from 
the XENONIO S2-only analysis [49]. We conservatively assume that the ionization yield 
drops to zero below 1.4 keVnr because the signal acceptance cannot be reliably estimated 
below this energy. The energy resolution is taken to be AEr = Er/E jiQy{Eji), while the 
acceptance is 95% on an exposure of 15 kg d. Yellin’s “Pmax” method [45] is used to set 
limits at 90% C.L., using the conservative assumption that all observed events are potential 
DM scatterings. Eollowing [43], we deduce Poisson upper limits from the results of the LUX 
experiment [6]. The background distributions for the LUX experiment have not been made 
public, so we cannot compute limits based on the profile likelihood ratio test statistic as 
done by the LUX collaboration. Because no events were observed in LUX in the low-energy 
nuclear recoil band, we instead use the 90% upper limit for an observation of 0 events for a 
Poisson process with no background, taking the efficiency from [6]. This simplification gives 
reasonable agreement with the limits reported by LUX in the low-mass region. 

In the left panel of Pig. 4 we compare limits from XENONIO, LUX, and SuperCDMS, 
for 10 GeV DM in the usual Umm-space plot. As pointed out in [43], the limits in Umm-space 
exhibit a simple scaling with the DM mass, 

v'^i^iEn) = ^v^UEr) ■ (4.1) 

Lnx' 

This can be visualized in a single plot, such as the right panel of Eig. 4, which shows the range 
of limits for SuperCDMS, when referred to 10 GeV DM using the scaling factor appropriate for 
Ge. Such a plot will contain, for every possible mass, the correct positioning of the CDMS-Si 
preferred region relative to SuperCDMS. Por sufficiently large DM masses, the exclusion limit 
from SuperCDMS drops below the lower boundary of the CDMS II Si contour, indicating 
strong tension between the inferred dark matter halos of the two experiments. The mass at 

computer code developed to determine these envelopes in general is available from the authors upon 
request. 

^^See [32, 35] for more details on reconstructing the published LUX results. 
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Vmin [km/s] [km/s] 

Figure 4. Exclusion limits from XENONIO (dot-dashed green), LUX (solid purple), SuperCDMS 
(dotted orange) overlaid on the preferred region (90% C.L.) for CDMS II Si (shaded blue) and best fit 
halo (dashed blue) in VminSp&ce. Left: Umin-space plot for 10 GeV DM. Right: range of SuperCDMS 
limits scaled with respect to 10 GeV DM using the factor appropriate for Si (shaded orange) (analogous 
to Fig. 7, right), overlaid with the preferred region for CDMS II Si (shaded blue). 


which the SuperCDMS exclusion drops below the lower limit of the CDMS II Si preferred 
region can be assessed using the methods of Sec. 3.2. 

The same information is more directly accessed in the corresponding momentum-space 
plots, shown in Fig. 5, which have the advantage of not requiring a choice of reference DM 
mass. Because the lower boundary of the CDMS II Si region and the SuperCDMS exclusion 
have some overlap along the /i-axis, their comparison corresponds to case (b) in Fig. 2. 
Using equation (3.14) and the minimum ratio between the SuperCDMS exclusion and the 
lower boundary of the CDMS II Si region, J = 0.22, we conclude that SuperCDMS excludes 
the CDMS II Si preferred region for all DM halos whenever rriy. >11 GeV. This corresponds 
to the dot-dashed red line in Fig. 3, left. 

Since Mq^ > Msi, we can applying the escape velocity constraints from Sec. 3.3. Using 
PB,max — 2.1 X 10“^ GeV/c, we find that SuperCDMS excludes CDMS II Si for all DM masses 
and all halos with escape velocity less than Vesc — 410 km/s. Indeed, this can be seen in the 
left panel of Fig. 4: for rriy. = 10 GeV, the lower boundary of the oreferred envelope contains 
points with Vmin > 410 km/s. Consequently, while the SuperCDMS exclusion does not 
cross the lower boundary of the envelope, a DM interpretation is still excluded for halos with 
sufficiently small escape velocities, since going to smaller DM masses would push the envelope 
further to the right. This implies that a consistent interpretation even for ^ H GeV 
will require scattering from DM particles in the tail of the DM velocity distribution above 
410 km/s. Of course, the Galactic halo almost certainly has a larger escape velocity than 
410 km/s, but we regard this example as an interesting proof-of-principle for a method which 
helps point towards particular DM mass ranges and halo properties. 
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Figure 5. Momentum-space plot comparisons between CDMS II Si and SuperCDMS. Left: exclusion 
limits from SuperCDMS (dotted orange) in p^j-space overlaid with preferred region of CDMS II Si 
(shaded blue) and best fit halo (dashed blue). Right: i?-ratio of SuperCDMS exclusion with lower 
boundary of CDMS II Si preferred region. Deviations from monotonicity are smoothed so that the 
limits h{pp) can be inverted to form R{h), as in equation (3.7). The minimum value of ^ = 0.22, 
implies that SuperCDMS excludes CDMS II Si across the full range of Vmin only for ^ II GeV. 
For halos with escape velocity less than Vesc — 410 km/s, SuperCDMS excludes CDMS II Si for all 
DM masses. The critical horizontal line at is above the top of the plot. 


5 Analysis Recipe 

In order to develop a halo-independent momentum-space analysis framework, we have had 
to introduce a good deal of notation, as well as several types of plots which were useful for 
pedagogical purposes. However, in practice the momentum-space method is straightforward, 
and here we summarize the approach. The minimal analysis recipe is as follows: 

1. For each experiment, compute limits (for null results) or best-fit regions (for signals) 
in pfi — h space as in Fig. 5, left panel. This is accomplished by constraining the 
differential event rate in Eq. (3.3) using whichever statistical technique is desired. For 
each individual experiment, this represents the universal halo-independent momentum- 
space result which can be derived from the data and this plot is the same for all DM 
masses, i.e. it is mass-independent.^^ 

2. To compare the consistency of two experiments, construct the R{h) plot from the ratio 
of the h{pR) curves for each experiment as in Fig. 5, right panel, using Eq. (3.7). The 
value of 6 (defined in Eq. (3.12)) extracted from such a plot determines the endpoint of 
the range of masses for which the two experiments are inconsistent, using Eq. (3.13). 

3. As an optional extension, if the target mass for the null experiment is greater than the 
target mass for the experiment which observes a signal, determine the value of Vesc for 

definition this contains all of the information required to make the standard Vmin-space plot for any 
specific choice of DM mass if desired. 
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which the two experiments are inconsistent for all halos and all masses, using Eq. (3.17). 
The value of PB,max required can be read off the momentum-space plot constructed in 
step 1. 

Note that in the recipe we have described it is never necessary to pick a fiducial and actu¬ 
ally construct the mapping into Umm-space. In particular, the plots in Fig. 4 are not necessary 
for the analysis, and were provided here only for pedagogical reasons. The momentum-space 
plot, which is unique and independent of the DM mass, contains all of the information 
necessary to construct the halo-independent VminSpace plot for any DM mass. Thus, in 
the presentation of new experimental data, a single momentum-space plot contains all of 
the information required for external groups to perform halo-independent analyses. In this 
way one automatically avoids introducing the errors in interpretation which may arise from 
those external to an experimental collaboration attempting to recreate quantitative details 
of experimental analyses. 

6 Conclusions 

In this paper we have shown how a simple change of variables, from Umm-space to recoil 
momentum space, allows a presentation of halo-independent DM direct detection results 
without a fiducial choice of the DM mass. While previous applications of halo-independent 
methods required a choice for m^, resulting in a proliferation of plots, the positive or null 
results of an experiment for all DM masses can be shown on a single pR — h plot. To 
compare results between experiments, an accompanying R{h) plot which shows the ratios of 
the various experimental constraints may also be useful. Such a plot is easily derived from 
the PR — h plot, although we emphasize that it this additional plot is not necessary for the 
analysis of a single experiment. Depending on the mass orderings of the detector materials, 
one can in some cases determine directly from single pR — h plot, or in all cases from the 
R{h) plot, whether a DM interpretation for any possible halo is ruled out either (a) for all 
DM masses, or (b) for a limited range of DM masses, where the upper or lower endpoint of 
this range can be computed from an equation like Eq. (3.14) and read directly off the R{h) 
plot. 

In our discussion, only certain simplified scenarios have been sketched, and more com¬ 
plex exclusion curves and preferred regions are possible in the momentum-space plot. How¬ 
ever, the generalization of our results is completely straightforward, and our sample analysis 
in Sec. 4 (and App. A) shows that our method can be easily applied to real data. As the 
momentum-space plot is unique, in that it does not need to be recreated for multiple DM 
mass values, we believe that it is ideally suited for any halo-independent analysis by an ex¬ 
perimental collaboration. This would allow the presentation of experimental results in a way 
which is as complementary as possible to the m — a plots already in common use. 
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Figure 6. Some possibilities for the variable Rab defined in Eq. (A.l). After rescaling to Vmin- 
space, a point satisfying Rab,x < 1 {Rab,x >1) signifies a stronger constraint from experiment A 
(B). This rescaling is given by F^, defined in Eq. (3.10), which satisfies the inequality Eq. (3.11). 
Thus any curve lying entirely below Rab = 1, such as curve (a), implies that for all DM masses the 
halo-independent Umm-space constraint will be stronger for experiment A than experiment B. The 
converse is true for any curve that lies entirely above Rab = Ma/Mb, such as curve (b). See text for 
a detailed discussion of the other curves and general cases. 
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A Comparisons Between Null Results 

While the ultimate purpose of the mass-independent method is as a diagnostic tool in com¬ 
paring emerging hints of DM detection with null direct detection results, in this appendix 
we consider the comparison between momentum-space exclusion curves for two different ex¬ 
periments. 

A.l Ratio Plots for Null Results 

We consider two direct detection experiments A and B, with target nuclei of mass Ma and 
Mb respectively. Since we are comparing two null results we can assume, without loss of 
generality, that Ma > Mb-^^ 

As discussed in Sec. 3.2, in some simple cases it is possible to determine the relative 
strengths of bounds from two different detectors directly from their h{pB,) plots. For instance, 
in Fig. 1 experiment A finds the strongest constraints in p^-space, and as any mapping to 
Umin-space will shift the curve for B further to the right relative to A, we may conclude that 
experiment A finds the strongest exclusion for all DM masses. However, in general the two 
exclusion curves may in principle cross any number of times while remaining consistent with 
the monotonically decreasing behavior. 

To treat the general case we again exploit the fact that mapping from momentum- 
space to Umin-space only shifts curves relative to each other in the horizontal direction. 
Each momentum-space exclusion curve is given by the function h{pB)] since this curve is 

^■^Note that this was not the case in comparing a null result to a putative positive result. 
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Figure 7. Exclusion limits from XENONIO (dot-dashed green), LUX (solid purple), and SuperCDMS 
(dotted orange) in ?;mm-space. Left: iimin-space exclusions for 10 GeV DM. Kinks and slight devi¬ 
ations from monotonicity are present in the SuperCDMS limit, occuring when the “Pmax” method 
switches the interval used to calculate the exclusion. Similar discrete artifacts are not present in LUX 
because an observation of zero events is used to calculate the limit, nor in XENONIO because the 
experiment observed a large number of events. Right: limits with Vmin calculated for a range of DM 
masses and scaled to = 10 GeV using the factor appropriate for Ge iiGe,xl^^Ge,x' i where /ice.x 
is the reduced mass for 10 GeV DM and Ge target. This has the effect of mapping SuperCDMS 
limits to a single line (dotted orange), independent of mass, while limits from Xe experiments are 
mapped to different locations depending on the DM mass. Dark bands correspond to the range of 
limits with masses 5-100 GeV, while the light bands show the range of limits for all possible masses, 
which facilitates comparison of the relative strengths of the limits for different DM masses. 


monotonically decreasing, we may again define the inverse mapping PR^h) as in Eq. (3.6). 
We now define the ratio of the momentum-space exclusions for the two experiments 


RAB{h) 


PR Ah) 
PR Ah) 


(A.l) 


which is the analogue of R{h) defined in Eq. (3.7). The subscript AB is meant to indicate 
the mass ordering, Ma > Mr. 

It is useful to consider first the general form of the Rab curves. In general, two ex¬ 
perimental exclusion lines in momentum-space may cross any number of times, and thus 
Rab may cross unity and change gradient any number of times. In particular, RabA need 
not be monotonic. Its asymptotic behavior is somewhat simpler, though. As —>■ oo, the 
asymptotic value of RabA is the ratio of the low-energy thresholds of experiments A and 
B. The behavior in the other extreme —)• 0 is determined by the constraint which extends 
to lower h at high pR] if A’s constraint extends to lower h, then RAB{h —)• 0) —)• 0, otherwise 
RAB{h —)• 0) —>■ oo. These two cases inform the definition of 6ab, the analogue of 6 defined 
in Eq. (3.12): 


r R'XAh) - 1, if RAB{h ^ 0) ^ oo 

1 R^fCh) - 1, if RAB{h ^ 0) ^ 0 
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Here, 5ab is always positive since we have assumed Ma > Mb- For curves that have 
RAB{h —)• 0) —>• 0, 6ab is the furthest distance the curve gets above 1, and for RAB{h -A 
0) —)■ oo, 5ab is the distance of closest approach to 1. 

With these general discussions in hand we now discuss the various classes of behaviour 
for the Rab curves, illustrated in Fig. 6. As in Sec. 3.2, the rescaling of Rab into Vmin-space 
is Rab,x = RabFx^ with = ^bx/ l-^Ax satisfying Mb/Ma < F^ <1. We first discuss the 
cases in Fig. 6, left: 

(a) Since Rab < 1 everywhere then for all DM masses A would give the strongest constraint 
over all of Umj,i-space, as B would shift further to the right than A in all cases. This is 
the case shown in Fig. 1. A gives the strongest constraint over all of Umm-space, and 
thus for all possible DM halos for all DM masses. 

(b) This is the complementary case to (a) as Rab > Ma/Mb everywhere. The rescaling 
in going to Umin-space is insufficient to alter the relative strength of the constraints, 
for all DM masses. In going to Umm-space, the smallest factor by which the curve (b) 
is rescaled is Mb/Ma, corresponding to M^ -A oo. This cannot push any point on (b) 
below unity, so in this case B gives the strongest constraint over all of Umm-space and 
hence for all possible DM halos and all DM masses. 

(c) This case corresponds to a situation where, for all DM masses, the two experiments 
provide complementary constraints. There is no DM for which either experiment will 
give the strongest constraints over all of Umm-space: at low h, A is a stronger bound, 
and at high h, B is a stronger bound. 

(d) The same as (c) but with the roles of A and B reversed: at low h, B is a stronger 
bound, and at high h, A is a stronger bound. 

Cases (b') and (c') in Fig. 6, right, are those for which A oi B can give the strongest 
constraint for a range of DM masses. As in Sec. 3.2, the endpoints of these ranges are 
determined by Mab{Sab), defined in Eq. (3.13). 

(b') If the DM mass mass satisfies 

M^ < Mab{Sab) , (A.3) 

B will give the strongest constraint over all of Umm-space and hence for all DM halos, 
(c') If the DM mass satisfies 

M^>Mab{Sab) , (A.4) 

then A will give the strongest constraint over all of Umm-space and hence for all DM 
halos. 

A.2 A Sample Analysis: Comparing LUX, XENONIO, and SnperCDMS 

The dark matter mass scaling relations introduced in section A.l for comparing null results 
can be used to concisely illustrate the mass ranges for which Xe-based and Ge-based exper¬ 
iments give stronger results. In the left panel of Fig. 7 we compare limits from XENONIO, 
LUX, and SuperCDMS, for a 10 GeV DM. As was demonstrated in Sec. 4, this can be visu¬ 
alized in a single plot, such as the right panel of Fig. 7, which shows the range of limits for 
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Figure 8. Comparison of mass-independent exclusion limits. Left: momentum space exclusion limits 
from XENONIO (dot-dashed green), LUX (solid purple), and SuperCDMS (dotted orange). Right: 
ratio of exclusion in p^-space of LUX to SuperCDMS (solid purple) and XENONIO to SuperCDMS 
(dot-dashed green). Deviations from monotonicity are smoothed so that the limits h^pu) can be 
inverted to form Rxe-Geih), as in equation (A.l). The minimum value of i5 = 0.25 implies that 
SuperCDMS is stronger than XENONIO across the full range of Vmin only for < 59 GeV. The 
SuperCDMS and LUX limits, on the other hand, probe complementary ranges of Vmin for all DM 
masses. 


each experiment, when referred to 10 GeV DM using the scaling factor appropriate for Ge. 
Such a plot will contain, for each possible mass, the correct positioning of LUX or XENONIO 
limits relative to SuperCDMS. Although one cannot directly infer the value of Vmin at which 
curves cross at each DM mass, it is manifestly clear from the plot that SuperCDMS is al¬ 
ways stronger than XENONIO at small DM masses, and that SuperCDMS and LUX are 
complementary across all DM masses due to the higher Vmin threshold of LUX. 

Fig. 8 shows the momentum-space exclusion plots for SuperCDMS, XENONIO, and 
LUX. The left panel shows the exclusions as a function of recoil momentum, while the right 
panel shows the ratio of LUX and XENONIO limits to the SuperCDMS limit, as a function of 
h. Comparing Fig. 8 to the cases enumerated in Fig. 6 indicates that the pair of limits from 
SuperCDMS and LUX corresponds to case (c), while the pair of SuperCDMS and XENONIO 
corresponds to case (b'). In case (c), both limits are complementary in Vmin space for all 
possible DM masses, and this is clearly apparent for SuperCDMS and LUX in Fig. 7. In case 
(b'), SuperCDMS is stronger than XENONIO for DM masses below the bound in equation 
(A.3). Using the value of h = 0.25 read off from the Rxe-Ge plot constructed from the 
SuperCDMS and XENONIO limits, we find that SuperCDMS is stronger for all Vmin when 
< 59 GeV. This behavior is also qualitatively visible in the right panel of Fig. 7. 
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